Complex magnetic behavior and high spin polarization in Fe3_2;Mna;Si alloys 
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FeaSi is a ferromagnetic material with possible applications in magnetic tunnel junctions. When 
doped with Mn, the material shows a complex magnetic behavior, as suggested by older experiments. 
We employed the Korringa-Kohn-Rostoker (KKR) Green function method within density-functional 
theory (DFT) in order to study the alloy Fea-aiMnaiSi, with < a; < 1. Chemical disorder is 
described within the coherent potential approximation (CPA). In agreement with experiment, we 
find that the Mn atoms align ferromagnetically to the Fe atoms, and that the magnetization and 
Curie temperature drop with increasing Mn-concentration x. The calculated spin polarization P at 
the Fermi level varies strongly with x, from P — —0.3 at x = (ordered FeaSi) through P = at 
X = 0.28, to P = +1 for X > 0.75; i.e., at high Mn concentrations the system is half-metallic. We 
discuss the origin of the trends of magnetic moments, exchange interactions. Curie temperature and 
the spin polarization. 

PACS numbers: 75.50.Bb, 71.20.Be, 71.70.Gm, 71.20.Lp 



I. INTRODUCTION 

Magnetic internietallic alloys show very rich physics 
depending on the degree of doping and chemical disor- 
der, which can therefore be used as "control parameters" 
allowing to tune the electronic and magnetic structure for 
desired effects. In particular, physical properties that are 
fundamental for technological applications in spintronics, 
such as the magnetization M, Curie temperature Tc, or 
spin polarization P at the Fermi energy, vary strongly 
with respect to these control parameters 

The alloy Fe3_a;Mna;Si, belonging to the wider class of 
Fea-i^TMjjSi alloys with TM a transition-metal elementfi 
is an example of such dependence on the degree of 
doping 1^ As the Mn concentration increases, the mag- 
netization of Fe3_2;Mn2;Si drops continuously from 5 fiB 
to zero; its temperature-dependent magnetic properties 
change from high-Tc (~800 K) ferromagnetism, through 
low-Tc ferromagnetism with re-entrant behavior at 70 K, 
to complex non-collinear magnetism; its calculated spin 
polarization increases from —30% to the ideal, half- 
metallic +100%^ and then drops again due to the non- 
collinear behavior. These observations are not new, 
however, there is a recent revival of the interest in 
Fe3_2:Mn^Si due to potential applications in magnetic 
tunnel junctionsi^ 

Motivated by this revival, we present here a theoret- 
ical study to the electronic and magnetic properties of 
Fea-j^Mn^^Si for < x < 1 based on ab initio calcula- 
tions. The choice of concentration range is motivated by 
the specific site preference of Mn for x < 1, so that the 
resulting state is ferromagnetic as we see later; for a; > 1, 
non-collinear magnetic configurations can occur. We pro- 
vide an interpretation of the magnetization drop as a 



function of concentration in terms of wavefunction sym- 
metry and hybridization, together with the requirement 
for local charge neutrality. We further propose that the 
increase in spin polarization up to the half-metallic point 
is due to the same mechanisms that cause the magneti- 
zation drop. Moreover, after extracting exchange inter- 
actions from the ab initio results, we calculate the Curie 
temperature using a Monte Carlo approach, and are able 
to reproduce the drop of Tc as a function of Mn con- 
centration. Finally, we discuss where our results do not 
agree with experiment, and we propose a possible rea- 
son for the disagreement; this is particularly the case for 
the re-entrant behavior and the value of magnetization 
at high Mn concentrations. 

The paper is structured as follows. In Sec. |n]we sum- 
marize the experimental and theoretical background on 
Fea-j^Mn^^Si. Sec. |in]is devoted to the description of our 
calculational approach. We continue with a presentation 
of our results on the magnetization and spin polariza- 
tion in Sees. IIVI and |V] and of the Curie temperature in 
Sec. IVIl In Sec. IVIII we discuss the limitations of our ap- 
proach, their consequences, and possibilities for a more 
accurate description. Our conclusions are summarized in 
Sec. lVml 



II. EXPERIMENTAL AND THEORETICAL 
BACKGROUND 

A. Experimental 

Quite a few experimental studies have been done on the 
magnetic properties of Fe3_a;Mna;Si, revealing a highly 
complex magnetic behavior dependent on the Mn con- 
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FIG. 1: Geometrical structure of FegSi and Fe3_a;Mna;Si. The 
A and C sites are occupied by Fe, while the B site is occupied 
either by Fe or by Mn. Square symbols show the positions 
of the Si atoms. In the left panel we show the octahedral 
coordination of the B-site to A and C Fe neighbors, while in 
the right panel we show the tetrahedral environment of the A 
and C sites. 



centration x. Here we recall the main results of these 
experiments, with emphasis on the concentration range 
< a; < 1 which interests us in the present work. At high 
concentrations, one reaches the Fe-doped MnsSi com- 
pound, which exhibits more complicated properties; e.g., 
MnaSi is an incommensurate antiferromagnet with a Neel 
temperature of about 25 K, while FesSi is a ferromagnet. 

Structure and site preference. FeaSi crystallizes 
in the DO3 structure consisting of a fee lattice with four 
basis atoms (see Fig. [1]). These are placed at (0,0, 0)a 
(Si atom), (i,i,i)a (FcA-atom), (i,i,i)a (Fee-atom), 
(|,|,|)a (Fee-atom), where a is the lattice constant. 
The Fca and Fee atoms are tetrahedrally coordinated to 
four Fcb and four Si atoms and exhibit equivalent elec- 
tronic ground-state properties due to symmetry. The Fee 
atoms are octahedrally coordinated to eight FeA,c atoms; 
their electronic properties therefore resemble somewhat 
bcc Fe, as we will discuss. 

When Mn is doped into FcsSi, it substitutes Fcb atoms, 
as is found by experiment. This appears to be part of 
a general trend found expcrimentallyii^i^ and modeled 
theoretically^ in which transition-element atoms which 
are to the left of Fe in the periodic table prefer to reside 
at the B-site when doped into FeaSi, while transition- 
element atoms which are to the right of Fe prefer to sub- 
stitute the Fe atoms at the A and C sites. 

Yoon and Booth^i^ report that, in the range < a; < 
0.75, the Mn atoms substitute the B-site Fe, as verified 
by hypcrfinc field measurements of Niculcscu et al.i^ 
Above X = 0.75, an increasing fraction of Mn atoms 
tends to reside at the A,C sites with equal probabil- 
ity. At a; = 1, this "swap" fraction is 12% according 
to Ziebeck and Webster— and 15% according to Yoon 
and Bootb^; i.e., even Fe2MnSi contains a small amount 
of chemical disorder. Fe3_a;Mnj;Si has thus the L2i, or 
fuU-Heuslcr, structure. Similar results are found in more 
recent hypcrfinc-field experiments^^ 



The measured lattice parameter changes linearly and 
only slightly as a function of the concentration x (see 
Ref i) from 5.653 A at a; = to 5.663 A at a: = 1, i.e., 
by about 0.2%. The linear change continues for higher 

X. 

Magnetic moments. As was shown by Booth et alt^ 
and Yoon and Booth,— the saturation magnetization M 
drops linearly from about 4.8 //b to 2.6 /iB per formula 
unit in the range < x < 0.75. The local magnetic 
moments depend strongly on the site. The B-site has a 
high moment of about 2.3±0.3 /iB between < a; < 0.75, 
which then drops gradually and vanishes at a; = 1.75; at 
the A,C sites (containing only Fe for low a;), the moment 
decreases from 1.4 /xb to 0.3 /iB as a; increases from to 
0.75. The net result is the aforementioned drop of the 
total magnetic moment. 

Anomalous temperature dependence of the 
magnetization. In the concentration range 0.75 < a- < 
1.75 a re-entrant behavior of the magnetization curve 
A/(T) is found>^i^ for temperatures from T = K up 
to the re-entry temperature Tr « 70 K, M{T) is in- 
creasing, while for T > Tr M[T) is decreasing, as ex- 
pected for a usual ferromagnet, up to the Curie tem- 
perature Tq. Thus, two values of the saturation mag- 
netization can be defined, one (Mq) as the actual mea- 
sured value M{T = 0), and another (Moxtr) as the ex- 
trapolated value of M{T) from data taken for T > Tr; 
evidently, Moxtr > Mq. For x ~ 0.75, in fact, where 
this effect just starts to appear, the two values are close 
(A/o « 0.95Moxtr), and apphcation of a magnetic field 
of 14 kOe can further saturate the sample so as to reach 
the value Moxtr- The anomalous temperature dependence 
of M was also found by Ziebeck and Webster— for the 
Fc2MnSi alloy, as well as by Nagano and co-workersH 
and by Ersez et alf^ who analyzed the effect via neutron 
scattering. However, in a recent work where the samples 
were grown by molecular beam epitaxy;^ this anomalous 
behavior was not observed (up to a; = 1.4). 

Curie temperature. The Curie temperature drops 
as a function of the Mn concentration a;, from about Tq sa 
800 K at a; = to Tc = K at a; « 1.75; i.e., for 
X > 1.75 no saturation magnetization is found, while the 
sample is at a complex non-collinear magnetic state. At 
X = 0.75, where the re-entrant magnetic behavior sets on, 
Tc = 375 K and Tr « 40 K, while at a- = 1, Tc « 220 K 
and Tr « 70 

Thermodynamic properties. The anomalous be- 
havior at Tr also shows in measurements of thermody- 
namic quantities. In particular. Smith et alJ^ have mea- 
sured a specific heat anomaly at Tr, most pronounced 
for the a; = 1 compounds. Furthermore, Miles et alt^ 
report a sharp peak of the thermal expansion coefficient 
at Tr = 60 K for the Fe2MnSi alloy {i.e., x = 1), while 
they find no such pronounced behavior at Tc. 

Heat treatment. As reported in Ref 0, the sam- 
ple preparation included a 24-hour heating at 830° C and 
water quenching; a different heat treatment (21 days at 
550°C and slow cooling) for samples with 0.95 < a; < 1.25 
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appeared to increase Mq, bringing it closer to Moxtr, 
but had no effect on the values of Tr or Tq. Ziebeck 
and Webster— and Smith et ali^ also used an annealing 
treatment at over 800° C for 24 hours water quenching. 
Miles et ali^ used two samples, one quenched from 800°C 
and one slowly cooled, with no change in the expansion 
coefficient and its anomaly at Tr. 



ranges between 10 and 20 mcV. Brown and co-workers2i 
analyze the behavior of the magnetic moments based on 
symmetry arguments and on a model by Swintendicb^ 
and conclude that, as the Mn concentration increases, 
the reduction of the exchange splitting leads to the drop 
of the FeA,c moment. 



B. Theoretical 

Only few theoretical results exist on the electronic 
structure of Fe3_j;Mnj;Si. Mohn and Supanetaii em- 
ployed an augmented spherical wave method and the 
local spin density approximation (LSDA) to density- 
functional theory to examine non-coUincar states in the 
ordered alloys Fes Si, Fe2MnSi (with the Mn atom at the 
B-site), FeMn2Si (with one Mn atom at the B-site and 
one at the A-site), and MnsSi. For FeaSi they found a 
ferromagnetic ground state, while non-collinear ground 
states were found for all other compounds. In partic- 
ular for Fe2MnSi, which is of interest here, Mohn and 
Supanetaii found a local energy minimum at the ferro- 
magnetic state, with a lower minimum for a spin-spiral 
along the [111] axis, at a q-vector of ^(0.45,0.45,0.45). 
According to their calculation, the ground state was an 
antiferromagnetic state, with the Mn moments alternat- 
ing along the [111] axis, while the moment direction was 
canted off the [111] direction by about 60°. The en- 
ergy difference between the non-collinear ground state 
and the ferromagnetic state in Fe2MnSi was reported to 
be around 0.8 mRyd (10.9 meV). A canting of the mag- 
netic moments below Tr was also assumed by Yoon and 
Booth^i^ in order to explain the neutron scattering data. 

More recently, ab initio calculations were presentecU^ 
on Fe3_2;Mn2:Si for < a; < 0.5 using the supercell 
method. Mn was considered to occupy the B site. It 
was found that the spin moments of Fe at the A and 
C sites are reduced in the presence of Mn nearest neigh- 
bors, which induces a redistribution of the FeA,c states; a 
drop of the total moment with increasing Mn concentra- 
tion was observed and attributed to the FeA,c moments. 
Furthermore, ab initio calculations on Fe3_a;Mna;Si and 
Fe3_2:MnSia; alloys were presented in Ref. [l^ the cal- 
culations here agree with the previous results that the 
magnetic moments of the A and C sites drop as a func- 
tion of Mn concentration, while it is found that the B-site 
atomic moments increase. The Fe moments appear to be 
higher than the Mn moments, so that the B-site average 
moment docs not change much. 

Further theoretical work appears in parallel with ex- 
periments. Szymahski et al^ examined the spin dynam- 
ics of Fe3_a;Mna;Si using neutron scattering at room tem- 
perature and at liquid nitrogen temperature, and fitted 
their results to effective interatomic exchange integrals 
which enter a Heisenberg Hamiltonian. The fitted values 
of the exchange constants depend on the number of neigh- 
bors considered; the nearest-neighbor exchange for Fc3Si 



III. METHOD OF CALCULATION, MODELS 
AND CUTOFF PARAMETERS 

Our electronic-structure calculations are based on 
density-functional theory within the generalized gradient 
approximation (GGA)^ to account for exchange and cor- 
relation effects. The local spin density approximation^^ 
was also used for comparison and proved to be inade- 
quate for the prediction of the correct magnetic ground 
state at low and intermediate Mn concentrations (see 
Sec. lIVp . Calculations were also performed within the 
"single-shot GGA";^ where, using the self-consistent 
LSDA spin density, Plsda, the total energy is calcu- 
lated within the GGA hmctional Eqqa[p]; i-e. one 
calculates i^cGA [plsda] • This approach is based on 
the idea that Plsda, as a trial density, is not too far 
from the self-consistent GGA density Pgga, so that, 
due to the variational character of the energy functional, 
£^gga[plsda] « -Egga[pgga]- The single-shot GGA is 
known, for instance, to correct the LSDA overbinding, 
giving an improved equilibrium lattice parameter, very 
close to the one predicted by the GGAi^ 

The Kohn-Sham equations are solved in most cases 
within the full-potential Korringa-Kohn-Rostoker (KKR) 
Green function method^ with exact treatment of the 
atomic cell shapes;^ using an angular momentum cut- 
off of Zmax = 3 and an integration mesh of 30 x 30 x 30 
points in the full Brillouin zone. The substitutional disor- 
der was described within the coherent potential approx- 
imation (CPA) . For the calculation of magnetically non- 
collinear states and static magnon spectra we employed 
the full-potential linearized augmented plane wave22r— 
(FLAPW) method as implemented in the FLEUR code^^ 
using a plane-wave cutoff of /cmax = 4 a.u.~^, an angular 
momentum cutoff of ^max — 8, mufhn-tin radii of 1.19 A 
for Fe and Mn and 1.222 A for Si, and an 17 x 17 x 17 
mesh in the full Brillouin zone. The FLAPW code was 
also used to cross-check the KKR results in some cases. 
Relativistic effects were taken into consideration within 
the scalar relativistic approximation, whereas spin-orbit 
coupling was not accounted for. 

Since the lattice parameter varies only slightly^ in the 
range < x < 1 we used the experimental value at a; = 0, 
a = 5.653 A, in all calculations. 

Magnetic excitations are modeled within a classical 
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Heisenbcrg Hamiltonian, 

H = -^Jijii ■ ij, (1) 

= -Y,.h,Nh-M, (2) 

where and Cj are unit vectors along the directions of 
the spin moments, Mi and Mj, at sites i and j, while 
the exchange pair-interaction constants Jij reflect the en- 
ergy cost for the mutual tilting of the moments. It is 
sometimes convenient for the discussion to use the form 
(21) with Jij — Jij/{MiMj). The constants Jij were ex- 
tracted from the spin-dependent KKR structural Green 
functions Gjj^-'(i?) and t-matrices tl''^\E) by virtue of 
the Liechtenstein formula^ 

= i-Im Tr p 4 (4 - tj) 4 (4 - tj) dE. (3) 

Here, Clj^'^^^E) and tl^^^E) are matrices in angular mo- 
mentum space and Tr denotes a trace over the angular 
momentum indices (Im). 

Having calculated the exchange constants, the Curie 
temperature of the compounds was calculated within a 
Monte Carlo method. For this purpose, exchange con- 
stants of atom pairs with distance up to 2.18 lattice 
constants were used; the simulation supercells included 
1536 magnetic atoms (512 unit cells). The Curie temper- 
ature was identified through the characteristic peak of the 
calculated susceptibility. The method, either in combi- 
nation with Monte Carlo simulations or with the random 
phase approximation for the solution of the Heisenberg 
model, has proven useful for the calculation of the ex- 
change constants and Curie temperature with a 10-15% 
accuracy in several cases, including elemental ferromag- 
nets and intermetallic alloys.—"— 



IV. GROUND STATE MAGNETIC MOMENTS 
AND CONFIGURATION 

The calculated magnetic moments (per atom and to- 
tal) as a function of concentration are depicted in Fig. [2] 
Here the Mn atoms were assumed to reside at the B site 
for all X. Evidently, the B-site atoms (Fe and Mn alike) 
are generally in a high-spin state, while the Fe atoms at 
the A and C sites arc in a low-spin state. A small, mono- 
tonic increase of the B-sitc atomic moment is observed 
as a function of the Mn concentration; the Fcb moment 
ranges between 2.6 and 2.9 fiB, while the Mn moment 
ranges between 2.1 and 2.6 /ie. In strong contrast, the 
FeA,c moments drop significantly as the Mn content in- 
creases, from 1.3 for pure Fes Si to 0.2 fiB for the 
ordered Fe2MnSi alloy. This drop of the FeA,c moment 
causes the decrease in the total magnetic moment per 
unit cell. 
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FIG. 2: (color online) Calculated magnetic moments per atom 
and total magnetic moment in the unit cell as a function of 
Mn concentration in Fe3_a;Mna:Si. 



This trend is in agreement with experimental 
findings^i^ii^ that the average moment at the B-site is 
high and remains more or less unaffected by Mn-doping, 
while the A and C-site moments drop significantly. (A de- 
viation from the experimental result is found for x > 0.8, 
where there is experimental evidence for re-entrant be- 
havior and reordering of spins.) The trend can be un- 
derstood by an analysis of the density of states and an 
understanding of the different wavefunction hybridiza- 
tion of MnB-FeA,c and Fcb-Fba.c atoms; we defer the 
discussion to section fVl 

In the calculation, the Mn dopants can be chosen to 
align ferromagnetically or antiferromagnetically to the 
FeaSi matrix. Total-energy calculations are then needed 
to identify the correct ground state, which experimentally 
is found to be ferromagnetic at least for concentrations 
X < 0.75. The calculations within the GGA show that 
the ferromagnetic state is stable for all Mn concentrations 
< a; < 1 (see Fig. [3|). However, the result within the 
LSDA is that the antifcrromagnctic state is more stable 
than the ferromagnetic one for x < 0.35, in clear disagree- 
ment with experimentr^ actually, for 0.2 ^ x < 0.8, the 
LSDA lowest energy is found in a disordered local mo- 
ment state of the form Fe3_xMn^_j,MnJSi, as the Mn 
atoms progressively change the moment orientation from 
antiferromagnetic (Mn-I-) to ferromagnetic (Mn'''). This 
indicates a possible non-collinear LSDA ground state 
(our calculations for the disordered alloys were always 
magnetically collinear). The dispute between the LSDA 
and GGA results at low concentrations was cross-checked 
and verified by a calculation within the FLAPW method, 
where a low Mn concentration of x ~ 0.125 was approxi- 
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FIG. 3: (color online) Results of total-energy calculations on 
the magnetic state of Fe3-2,Mna;Si. The energy difference 
AE = E{AF) — E(FM) between the antiferromagnetic and 
the ferromagnetic alignment of Mn atoms with respect to the 
Fee atoms is shown {AE represents here values per unit cell). 
The GGA predicts the correct ferromagnetic state, while the 
LSDA yields an antiferromagnetic orientation of the Mn mo- 
ments with respect to the Fes Si matrix at low concentrations, 
which turns to ferromagnetic at high concentrations. A single- 
shot GGA calculation improves the LSDA result somewhat, 
but not quite. 



mated by construction of a large supercell (the same lat- 
tice parameter was used in LSDA and GGA calculations, 
see Sec. [ml). 

The failure of the LSDA to predict the correct mag- 
netic ground state can be attributed to the exchange 
and correlation part of the total energy, rather than the 
single-particle energies. We arrive at this conclusion for 
two reasons. Firstly, the GGA density of states is very 
similar to the LSDA density of states (when both are 
calculated in the ferromagnetic configuration). Secondly, 
we attested our suggestion by using the single-shot GGA 
(described in Sec.lIIip. which changes only the exchange- 
correlation part of the total energy, while retaining the 
LSDA single-particle energies. The total-energy results 
are then improved significantly, although not entirely, to- 
wards the correct magnetic ground state as can be seen in 
Fig. 131 We note that there is no general rule favoring the 
GGA over the LSDA as far as the magnetic properties 
are concerned. For example, long-wavelength spin-wave 
spectra calculated within the adiabatic approximation 
agree rather well with experiment if the LSDA is used 
(see, e.g., Pajda et al^ for Fe, Co and Ni, or Buczek et 
ai— for intermetallic alloys). There are even reported 
cases, such as FeaAl?^ where GGA gives the wrong crys- 
tal structure, while LSDA corrects the structure as well 
as the magnetic moments. 



V. DENSITY OF STATES, MAGNETIC 
MOMENTS, AND HALF-METALLIC BEHAVIOR 

A. Trends of DOS with concentration; magnetic 
moments 

The density of states (DOS) of Fea-a^MuajSi at various 
concentrations x is depicted in Fig. m The gross features 
have been analyzed in the past in studies of full Heusler 
alloysi^ Most important points specifically for our dis- 
cussion are: (i) The hybridization of c?-states of t2g char- 
acter of the B-site atom (Fe or Mn) with t2g states of the 
A and C-site Fe atoms, (ii) The strong-ferromagnet char- 
acter of the B-atoms as opposed to the weak-ferromagnet 
character of the (A,C) atoms, (iii) The progressive shift, 
for charge-neutrality reasons, of the majority-spin states 
at the B-site as the Mn concentration is increased, drag- 
ging with them the (A,C)-site t2g states due to hybridiza- 
tion and affecting the Fca.c moment and DOS. We now 
discuss these points in detail; the trends of the magneti- 
zation and the appearance of half-metallicity are directly 
connected to this behavior and also discussed in this sec- 
tion. 

(i) The d-states of Fee and Mn are split by the octa- 
hedral environment into two irreducible representations: 
the t2g, including the dxy, dyz, and dxz orbitals, and the 
eg, including the dj.^_y2 and orbitals {x, y, and z are 
implied to be directions along the cubic crystal axes). 
The FeA,c rf-states are also split by the tetrahedral en- 
vironment. However, it is also important that d orbitals 
of FeA,c atoms interact with each other, although the 
Fca.c atoms are second-nearest neighbors. This interac- 
tion produces bonding-antibonding splittings between t2g 
and ti„-type of states, and also between Cg and e„-type 
of statesjS^ 

Hybridization between c?-states of B-site atoms with 
Fca.c atoms is allowed only among the states of t2g or 
the ones of Cg character; by symmetry, the tiu and e„ 
states remain oblivious to the d-states of their B-site 
neighborsi^ The symmetry-decomposed DOS is shown 
in Fig. [51 for FeaSi, and is completely analogous in 
Fcs-a-Mn^rSi, with the peaks appropriately shifted as dis- 
cussed below. 

(ii) On the one hand, the local DOS of the B-site atoms 
(either Fe or Mn) has in all cases the characteristics of the 
DOS of a strong ferromagnet. Namely, the B-site local 
DOS of one spin direction (here spin down) is very low at 
Ey and in a region around E-p , as is evident by inspection 
of Fig. m Consequently, on-site transfer between spin-up 
and spin-down charge is energetically expensive, because 
a strong shift of d-bands is involved. This, together with 
the requirement of local charge neutrality in a metal- 
lic system, stabilizes the local atomic moment of Fe or 
Mn at the B-site against perturbations (such as change 
of concentration or lattice parameter). Therefore there 
is only weak dependence of the Fe and Mn moment on 
concentration, as seen in Fig. [21 The strong-ferromagnet 
behavior of the B-site is favored by its octahedral bcc-like 
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FIG. 4: (color online) Spin- and atom-resolved density of states of Fes-^^MnajSi alloys for Mn concentrations < a; < 1. In all 
cases, the Mn atom taken to reside at the B-site. Up and down arrows indicate the majority and minority-spin, respectively. 
In (a) and (b), arrows indicate the positions of the DOS peaks near Ef which are responsible for the moment trends, as 
discussed in the text. The calculations were done within the KKR-CPA. In each plot, the upper panel corresponds to spin-up 
(majority-spin) DOS and the lower panel (with inverted ordinate) to spin-down (minority-spin) DOS. The orange-shaded area 
under the red line corresponds to Mn, the full blue line to Fee and the double green line to FeA.c- 



coordination, which results in large bonding-antibonding 
splittings via ^29-^23 and Cg-eg-hybridization of the B-site 
atoms with Fca.c neighbors. The line-shape of the local 
DOS at the B-site is reminiscent of bee iron. 

On the other hand, the local DOS of the FeA,c atoms 
has the characteristics of a weak ferromagnet, i.e., d- 
states of both spin directions are present at and around 
Ep. This allows for energetically cheap transfer between 
spin-up density and spin-down density, and therefore rel- 
atively easy change of moment. Note from Fig.[5]that the 
drop of the FeA.c-moment is strongest for 0.2 < a; < 0.4, 
at the same concentration range when the FeA,C"DOS at 
E-p is large for both spins in Fig. SI i.e., when the weak- 
ferromagnet character is most evident; before x = 0.2 
there is a DOS valley at Ep for spin-up, and after x — 0.4 
there is a valley for spin-down. The weak-ferromagnet be- 
havior of Fca.c is favored by the ti„ and states forming 
the spin-down peak of FeA,c around Ep in FeaSi. 

(iii) Now we consider the consequences of observations 
(i) and (ii) (see Figure El). We first focus on the B-site 
in Fea-j^Mnj^Si. Since a metal must show approximate 
local charge neutrality, the Mn states must be appropri- 
ately shifted with respect to those of Fe so that one less 



electron is accommodated by Mn. This is achieved by a 
shift of the spin-up t2g and Cg peaks, which are just under 
Ep for the Fcb atom, so that they fall at Ep for the Mn 
atom (as indicated by arrows in Fig. |3l^a,b) ) . Now, there 
is a peak of t2g character (also indicated by arrows in 
Fig. |4lja,b)) in the Fca.c DOS, associated by hybridiza- 
tion to the aforementioned B-site peak. This, at low Mn 
concentrations, is associated more to the Fcb DOS, while 
at high Mn concentrations it is associated more to the Mn 
DOS. At intermediate concentrations, as the Mn content 
increases, the FeA,c peak is dragged to higher energies 
and starts crossing the Fermi level, depriving the Fca.c 
atoms of spin-up charge. This is readily compensated (to 
maintain local charge neutrality) by a shift of spin-down 
FeA,c tiu states from Ep to slightly below Ep, gaining 
spin-down charge. The net effect is a reduction of the 
Fca.c spin moment, accompanied by the appearance of a 
spin-down gap at E-p, i.e., the signature of a half- metallic 
behavior. 

Thus we see that the ^29-^23 hybridization between 
states at the B and A,C sites, together with the require- 
ment of local charge neutrality, leads to the drop of Fca.c 
moment as the Mn concentration increases. 
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FIG. 5: (color online) Symmetry-resolved local density of 
states of FeaSi at the Fee (top) and Fba.c (bottom) atoms. 
In each plot, the upper panel corresponds to spin- up and the 
lower panel to spin-down DOS. 



In spite of the calculated drop of the magnetization per 
formula unit, the values are still too high compared to 
experiment,— although they are in agreement with pre- 
vious calculationsi^ii^'i^i^ For Fe2MnSi, the calculated 
value is 3 /ib, while the experiment gives Mq ~ 1.5 fiB 
and Moxtr ~ 2.3/iB. We propose an explanation of this 
discrepancy in Sec. IVIII 
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Spin polarization and transition to half-metallic 
behavior 



FIG. 6: (color online) Spin- and atom-resolved density of 
states of Fes Si (top) and Fe2MnSi (bottom) around E-e- Ar- 
rows indicate the shift of the peaks upon substitution of Fe 
by Mn that lead to the reduction of the moment and half- 
metallic behavior. 1. The B-site spin-up peak shifts higher in 
order to dispose of one electron. 2. The spin-up FeA,c peaks 
are dragged along due to hybridization, disposing of spin-up 
charge also at the FeA,c atoms. 3. The spin-down peak at 
FeA,c must shift lower in order to maintain local charge neu- 
trality. The result is a reduction of the magnetic moment and 
an increase of the spin polarization at Sp. In each plot, the 
upper panel corresponds to spin-up and the lower panel to 
spin-down DOS. The orange-shaded area under the red line 
corresponds to Mn, the full blue line to Fcb and the double 
green line to Fca.c- 



The spin polarization P at Ep is defined as 

p ^ n^iEp) - %(-Ef) 

with n^(i?F) and ni{Ep), respectively, the spin- up and 
spin-down DOS at E^- The mechanism described in the 
previous subsection, involving a shift of the peaks around 
E-p, leads to a drastic change of P as a function of the 
Mn concentration. This is visualized in Fig. [T] At low 
X, P is negative, approximately —0.3 (the negative sign 
means that the DOS at Ep is dominated by minority-spin 
carriers). At higher concentrations P crosses zero and 
reaches high positive values, as the spin-up DOS peak 
shifts to Ep, while the spin-down peak is retracted below 
Ep. At around x = 0.75 the spin polarization reaches the 
highest possible value of P = 1, and the alloy becomes 
half-metallic. At this point the (experimentally found) 
Curie temperature is about 370 K. A half-metallic be- 
havior of ordered Fe2MnSi was also found in previous 
calculations ^li^ii^i^^ 



This is a rare occasion in which a continuous change 
of a material parameter — here the Mn concentration — 
results in a continuous change of the spin polarization 
over such an extended range. Assuming that this effect 
is present in experiment, it could be efficiently used to 
compare to each other various experimental methods of 
probing the spin polarization (such as spin-polarized pho- 
toemission spectroscopy, positron annihilation, Andreev 
reflection, or tunneling magnetoresistance) . It should be 
noted that the half-metallic behavior is already present at 
X = 0.75, i.e., before the start of the re-entrant behavior. 



VI. EXCHANGE INTERACTIONS AND CURIE 
TEMPERATURE 

The exchange pair-interaction constants were cal- 
culated for several Mn concentrations as described in 
Sec. mil They are plotted as a function of distance in 
Fig. [5] In all cases, we find that the dominant contri- 
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FIG. 7: Calculated spin polarization P at Ef as a function of 
Mn concentration x in Fe3_a;Mna;Si. For x > 0.7 the alloy is 
half- metallic. The polarization of bcc Fe (dotted line) is also 
shown for comparison. 



bution comes from the first-neighbor interaction Ji be- 
tween the sitc-B atom (Mn or Fe) and FeA.c- The Fcb- 
FeA.c interaction, Ji(Fe-Fe), is between 20 and 10 mcV, 
depending on concentration, and the Mne-FcA.c interac- 
tion, Ji(Mn-Fe). is between 3-5 meV; next-nearest neigh- 
bor interactions (Fe-Fe, Mn-Mn or Fe-Mn) are typicahy 
at least one order of magnitude weaker. Therefore we 
expect the Curie temperature trend with concentration 
to follow the behavior of the averaged Ji , at least quali- 
tatively. 

The trends of Ji with concentration are depicted 
in Fig. [21 together with the moment of the Fba.c 
atoms, M(FeA,c)- There is a clear correlation between 
•/i(FeB-FeA,c) and M(FeA,c), while Ji(Mn-FeA,c) seems 
unaffected by the drop of M(FeA,c)- We now discuss 
these observations. The pair-interaction energies Eij ~ 
—Jij ii ■ ij, determined by the electronic structure, con- 
tain the absolute value of the atomic spin moments in a 
non-trivial way. By this we mean that, if the moments 
are varied by some external parameter {e.g., here, by 
changing the concentration), Eij can be affected either 
just by the variation of the absolute value of the mo- 
ment (as suggested by the alternative form derived from 
Eq. Eij = — Jij Mi ■ Mj), or also by an alteration of 
the exchange mechanism, which is induced by the change 
of the electronic structure via the external parameter and 
affects the constants Jij. 

Apparently, in Fe3_a;Mna;Si we are faced with both 
situations. On the one hand, the dominant trend for 
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FIG. 8; (color online) Pair exchange constants as a function 
of the inter-atomic distance. Left: 10% Mn concentration. 
Right: 70% Mn concentration. Note the different energy scale 
in the upper and lower panels. The lines are guides to the eye. 
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FIG. 9: (color online) Nearest neighbor exchange interactions 
Ji between Fee and FeA,c atoms and between Mna and FeA.c 
atoms, as calculated by the Liechtenstein formula ((3]). The 
FeA,c moments, M(FeA,c), are also shown, to demonstrate 
the correlation between the decrease of Ji and Af(FeA,c). 



Ji (Fe-Fe) comes from the reduction of the Fca.c mo- 
ment as the Mn concentration x increases, although there 
seems to be also an alteration of Ji, since the drop of 
M(FeA,c) is faster than the drop of Ji. On the other 
hand, Ji(Mn-Fe) is left practically unaltered despite the 
strong reduction of Af(FeA,c)- In order to interpret this, 
we again observe the shifting of the peak, indicated by 
an arrow in the Fca.c DOS in Fig. [H as a function of 
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X. This hybridization-induced peak coincides more and 
more with its associated Mn peak as the Mn concen- 
tration is increased. Since the Fermi level bisects the 
Mn peak (and increasingly more the Fca.c peak), the 
double- exchange mechanism sets in progressively more 
and more, favoring ferromagnetic alignment of the mo- 
ments. (We remind the reader that the double-exchange 
mechanism is present when half-filled states of the same 
spin hybridize with each other, resulting in a band broad- 
ening and a gain in energy; in a tight-binding picture, 
the kinetic energy is lowered by the inter-atomic hopping 
of electrons, allowed by the half-filled band.) Thus, the 
progressive shift of the indicated peak in the FeA,c DOS 
causes two competing effects: a reduction of the moment 
A'/(FeA,c) (as discussed in Sec. |V| and a strengthening 
of the Mn-Fe pair exchange interaction. These effects by 
and large cancel each other in the Heisenberg energy ex- 
pression, and the net result is only a weak dependence of 
Ji (Mn-Fe) on concentration. 

We close the discussion on the exchange parameters 
with the following comments on the calculations. In the 
present work, the exchange constants were calculated 
starting from the ground-state (ferromagnetic) config- 
uration. As a test, however, we calculated the Mn-Fe 
exchange constants starting from the antiferromagnetic 
(AF) state (i.e., with the Mn moments antiferromagneti- 
cally aligned to the Fes Si matrix) for a Mn concentration 
of X = 0.1, and compared with the result starting from 
the ferromagnetic (FM) ground state. As can be seen in 
Fig. [TOl there is a strong qualitative difference in the cal- 
culated value of the nearest-neighbor Mn-Fe interaction 
in the two cases: starting from FM, we obtain a ten- 
dency to retain ferromagnetism (positive exchange con- 
stant); while starting from AF, we obtain a tendency to 
retain antiferromagnetism (negative exchange constant). 
For the more distant neighbors the two calculations give 
quantitative, but not so much qualitative, differences. 
This discrepancy demonstrates the significant change in 
electronic structure at high angles, far beyond the as- 
sumptions of a Heisenberg model. The discrepancy is 
not observed at higher Mn concentrations, when the elec- 
tronic structure of the Fes-atoms at i^p is dictated more 
and more by the hybridization of their d-states with the 
Mn d-states. In addition, we note that, in the case of 
Fe2MnSi, the LSDA result for the Mu-Fca.c nearest- 
neighbor exchange parameter is weaker by approximately 
1/3 compared to the GGA result. This reflects also on 
the magnon spectra that are discussed in Sec. I VIII Such 
a discrepancy between LSDA and GGA does not show 
up for FcsSi. 

Calculated Curie temperatures Tc of Fcs-ajMn^jSi are 
shown in Table U together with experimental results. 
The experimental flnding of a reduction of Tc with in- 
creasing concentration is reproduced, although the cal- 
culated results systematically underestimate the exper- 
imental values. The reason for this trend is obviously 
the reduction of Ji(Fe-Fe) as a function of concentra- 
tion, in conjunction with the comparatively low values of 
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FIG. 10: (color online) Calculated Mn-Fe exchange interac- 
tion Jij via the Liechtenstein formula ((3| for a Mn concen- 
tration of a:: = 0.1, starting from two different configurations. 
Left: ferromagnetic starting point (i.e., with the Mn moments 
ferromagnetically aligned to the Pea Si matrix), which is also 
the ground state. Right: antiferromagnetic starting point 
(i.e., with the Mn moments antiferromagnetically aligned to 
the FesSi matrix). There is a strong qualitative difference in 
the dominant nearest-neighbor interaction, emphasized by an 
exclamation mark, while the more distant interactions change 
quantitatively, but not so much qualitatively. This discrep- 
ancy is not observed at higher Mn concentrations, when the 
electronic structure of the FeA.c-atoms at Ey is dictated more 
by the hybridization of their d-states with the Mn d-states. 
The lines are guides to the eye. 



TABLE I: Calculated Curie temperature of Fe3-2,Mni,Si at 
X = 0, 0.5, and 1. Experimental results (Ref. 0) are also 
shown. 



X 


Tc [K] (exp.) 


Tc [K] (calculated) 


0.0 


803 


730 


0.5 


450 


320 


1.0 


200 


160 (ordered; Mn at B) 






200 (12% Mn at A,C) 



Ji (Mn-Fe) which become important at high concentra- 
tions. Note that, at x = 1, the Fe2MnSi alloy is found 
expcrinicntallyiS to possess a degree of disorder in the 
form of a Mne-FeA.c swap of 12%. When considering 
such a swap in the calculations, we found an increased 
Tq of 200 K, mainly because the Fe atoms replacing Mn 
at the B-site have a stronger exchange interaction with 
the FeA,c neighbors. However, a possible Mn clustering 
at this concentrations, which could affect the value of 
the exchange interactions, cannot be taken into account 
within the CPA. This is discussed in more detail in the 
next section. 

Before closing this section, we show an interesting be- 
havior of the Mn sublattice magnetization at low or in- 



10 




FIG. 11: (color online) Monte Carlo simulation of the magne- 
tization and susceptibility x vs. temperature for Fe2.5Mno.5Si. 
It is evident that the Mn sublattice magnetization drops much 
faster than the Fe magnetization. The susceptibility peak sig- 
nals the Curie temperature. 



FIG. 12: (color online) Calculated flat spin-spiral disper- 
sion relations of Fe2MnSi along the [111] direction within the 
LSDA and GGA. The abscissa q corresponds to the point 
[q,q,q)2n /a of the Brillouin zone. The FLAPW method was 
used in this calculation. 



termediate concentrations. Kcpa et al^ and Ersetz et 
al^ find experimentally that the Mne moment at room 
temperature is significantly lower than what is reported 
close to T = K. E.g., for x = 0.55^ M(Mn) = 1.42 /^b 
at room temperature instead of 2.2 ^b- Therefore they 
propose that the Mn average moment drops with tem- 
perature faster than the Fe moment. This is reproduced 
by our Monte Carlo simulations, and is due to the weak 
coupling of Mn to Fca.c compared to the coupling of 
Fcb to FeA,c- As we see in Fig. [TTJ the Mn magnetiza- 
tion curve does not follow the critical-behavior form with 
an inflection point of the Fe curve, but rather drops al- 
most linearly. This is also reflected in the Mn sublattice 
susceptibility (not shown here), which does not show a 
peak ai T = Tq. 



VII. DISCUSSION ON THE RE-ENTRANT 
BEHAVIOR; LIMITATIONS OF PRESENT 
CALCULATIONS 

In Sec. In] we summarized what is known on the anoma- 
lous behavior of the magnetization of Fc3_^Mna;Si for x > 
0.75. The magnetization M (T) increases for < T < Tr 
(contrary to the behavior in a usual ferromagnet), then 
decreases again up to Tq; Tr varies with concentration, 
starting from low values at x = 0.75 and saturating at 
about 70-80 K a,t x = 1. This so-called re-entrant behav- 
ior and the resulting re-ordered phase, which has been 
the focus of many experimental works, could not be re- 
produced by calculations within the approximations used 
in the present work. 

There is experimental evidence^iiS based on neutron 
scattering data, that the re-entrant behavior arises from 
transverse ordering of the Mub magnetic moments. This 



means that, in the ground state, the Mn moments are 
partly canted to the direction of average magnetization. 
As the temperature is increased, the Mn moments ab- 
sorb energy by aligning their spins and thus the average 
magnetization increases. 

Mohn and Supanetaii performed first-principles 
(LSDA-based) calculations of non-collinear magnetic 
structures for spin-spirals of several wavevectors q. For 
ordered Fe2MnSi they found that the ferromagnetic 
state represents a local minimum, while a spin spiral of 
wavevector q = (27r/a)(0. 45, 0.45, 0.45) has a lower en- 
ergy by about 0.75 mRyd (10.2 meV). The absolute en- 
ergy minimum which they found is very slightly lower, for 
a state which shows a canted antiferromagnetic ordering 
of the Mn moments along the [111] axis, with a canting 
angle of 60° with respect to the [111] axis. 

We also performed calculations of magnetically non- 
collinear spin spirals along the [111] axis u sing the 
FLAPW method.^ However, the results of Ref. [I3 were 
not reproduced. Our results are shown in Fig. [T^ Cal- 
culations within both the LSDA and the GGA show 
a ferromagnetic ground state {i.e., a global minimum 
at (f = 0), while a local minimum is found close to 
q = (27r/a)(0. 45, 0.45, 0.45), i.e., at the point where the 
spin-spiral minimum is found in Ref. [TtI . The energy dif- 
ference between the two (local and global) minima is of 
the order of 11 mcV within the LSDA and 13 meV within 
the GGA. We do not know the origin of the discrepancy 
between our calculations and the calculations of Ref. Il3. 
A possible source of discrepancy is the use of a spherical 
potential approximation in Ref. ITtI as opposed to a full 
potential calculation here. 

However, these non-collinear calculations (ours as well 
as the ones of Ref. [TtI ) neglect a 12% Mn-Fe swap which 
is seen experimentally in Fe2MnSi; i.e., experimentally 
this alloy has a small degree of disorder. Given the small 
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calculated energy differences between the local and global 
minima (about 10 meV per formula unit), a correct de- 
scription of this swap can have important consequences. 
We attempted to model the swap by calculating the elec- 
tronic structure of (Feo.94Mno.o6)A,c(Feo.i2Mno.88)BSi 
within the CPA; here, the indexes (A,B,C) refer to the 
corresponding positions in the unit cell. The magnetic 
structure was subsequently investigated both by calculat- 
ing the exchange constants and the resulting configura- 
tion at T = K and by calculating a possible (collinear) 
disordered local moment state at the central site. How- 
ever, in all cases the outcome was a ferromagnetic ground 
state. 

What is missing within the CPA description is the 
short-range configurational order, i.e., the possibility to 
describe clustering of MnA,c atoms around Mub atoms. 
We speculate that this swap and clustering causes the 
canting of some Mne moments by interaction with the 
neighboring Mua.c atoms and is therefore essential for 
the appearance of the re-entrant behavior (such a sce- 
nario was already suggested in Ref . |^ . Possibly, as the 
temperature is increased, the canting of the Mub mo- 
ments is reduced on the average and the total magneti- 
zation increases; this hypothesis would require a weak 
coupling of the canted Mne moments to their Mua.c 
neighbors compared to the coupling of the non-canted 
Mub moments to their FeA.c neighbors. 

Our hypothesis is supported by the following exper- 
imental findings, (i) The re-entrant behavior and the 
swap appear in the same concentration range (x > 0.75). 
(ii) The re-entrant behavior is sensitive to the heat treat- 
ment of the alloy; after proper annealing, it was found 
that the re-entrant behavior smoothens, although Tr 
does not changei^ This, in conjunction with the calcu- 
lations showing that the un-swapped state is the ground 
state, suggests that annealing causes a fraction of the Mn 
atoms to return from the A,C site to the B-site, so that 
the number of Mn clusters (and canted Mn moments, as 
we suspect) lessens. If this hypothesis is true, Tr should 
be indeed unaffected by annealing, because it would cor- 
respond to a finite-size effect (characteristic exchange- 
energy scale of a small cluster) rather than a phase tran- 
sition. Such a possibility has been suggested by Nielsen 
and collaborators?^ based on Monte Carlo calculations 
of model systems, (iii) Clustering would result in a local 
environment which is closer to MuaSi, which is known to 
show antiferromagnetic behavior. 

Furthermore, it is observed that the "smoothening" of 
the re-entrant behavior after annealing is accompanied 
by an increase of both Mq and Afextr (shown in Fig. 7 of 
Ref. [3 for Fe1.75Mn1.25Si). We therefore make the plausi- 
ble assumption that this behavior is present also at con- 
centrations X < 0.75 < 1, and that the ground state, 
with all Mn atoms being at the B-sitc, will show a higher 



magnetic moment than Mextr- This could resolve the dis- 
crepancy between the measured^ magnetization value of 
A/oxtr ~ 2.3 /iB compared to calculated value of 3 /iB in 
FeaMnSi. 

First-principles investigations of such short-range- 
order effects require an approach beyond the CPA, 
e.g. by use of the non-local CPA^ or large-supercell 
techniques^ and we therefore defer this study to a fu- 
ture work. 



VIII. CONCLUSION 

The electronic and magnetic structure of the magnetic 
intermetallic alloy Fe3_a;Mn2:Si, for < a; < 1 has been 
investigated using density-functional theory within the 
GGA, together with the CPA to describe disorder. We 
find that important experimental findings, such as the 
trends of the magnetization and of the Curie temperature 
as a function of concentration, are reproduced. They can 
be interpreted in terms of single-particle energies with 
the help of the density of states, by using simple physical 
arguments, namely symmetry-dependent hybridization 
of wavefunctions and local charge neutrality. Quantita- 
tively, the drop of the magnetization as a function of Mn 
concentration is underestimated, especially at a; > 0.75, 
however this could be due to Mn-Fe swap and Mn cluster- 
ing in experiment that cannot be captured by the CPA; 
the same applies to the re-entrant behavior. It is there- 
fore worthwhile to investigate Fea-xMn^jSi beyond this 
approximation in the future. 

However, it is rather surprising that the two most com- 
mon approximations to DFT, namely the LSDA and the 
GGA, give ground states that differ at low concentration 
not only quantitatively but also qualitatively, although 
they almost agree on the single-particle spectrum. The 
superiority of GGA that we find here is common but cer- 
tainly not general, which shows the need for theory to go 
hand-in-hand with experiment for the understanding of 
magnetic intermetallic compounds. 

Concerning the relevance to spintronics, the most im- 
portant finding here is the continuous variation of the 
spin polarization at E-p over a wide range, between 
—0.3 <P<lforO<a:;< 0.75, i.e., in the region where 
theory and experiment are in reasonable agreement and 
before the onset of the re-entrant behavior. Since the 
polarization at E-p is a property that is notoriously dif- 
ficult to measure with precision, the variation which is 
found theoretically could be used to improve or calibrate 
the methods of measurement of P in a single type of 
material, so that spurious effects in measurement can be 
treated on the same footing and understood better. 
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